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We discuss the effect of quantum interference on transport through a quantum dot system. We 
£NJ . introduce an indirect coherent coupling parameter a, which provides constructive/destructive in- 

terference in the transport current depending on its phase and the magnetic flux. We estimate 
the current through the quantum dot system using the non-equilibrium Green's function method 
as well as the master equation method in the sequential tunneling regime. The visibility of the 
Aharonov-Bohm oscillation is evaluated. For a large inter-dot Coulomb interaction, the current is 
strongly suppressed by the quantum interference effect, while the current is restored by applying an 
oscillating resonance field with the frequency of twice the inter-dot tunneling energy. 

PACS numbers: 

'd 

1 I. INTRODUCTION 

Quantum phase coherence in mesoscopic systems is strikingly demonstrated with the principle of superposition, or 
interference experiments. Aharonov-Bohm (AB) interference is the most fundamental type and has been experimen- 
tally confirmed in metallic and semiconductor rings. Recently, in interference experiments with an Aharonov-Bohm 
ring containing a quantum dot (QD) in one of the arms, quasi-periodic modulation of the tunneling current has 
' been demonstrated as a function of the magnetic flux through the ringi^. This confirms that phase coherence is 
maintained during the tunneling process through a QD. The Fano effect is another type of interference in mesoscopic 
Q ■ physics, which occurs in a system in which discrete and continuum energy states coexist^. 

More recently the AB oscillations of a tunneling current passing through a laterally coupled double quantum dot 
(DQD) system were observed^. These experimental results have motivated theoretical investigations of electron 
transport through such a systen&2ii£. DQD has been attracting attention as an important device structure for 
entangled spin qubit operationsiiii2*A£. There is also an interesting theoretical prediction that cotunneling currents 
passing through spin-singlet and triplet states have different AB oscillation phases^. 

In this paper, we consider the transport through an AB interferometer containing a laterally coupled DQD. We 
introduce the indirect coupling parameter a, which characterizes the strength of the coupling via the reservoirs between 
two QDsi£. A system with the maximum coupling \a\ — 1 has already been widely studied theoretically^*^. In 
actual systems, however, such a case is very special and most experimental situations correspond to \a\ < 1. The 
situation where a = has also been explored in the context of the orbital Kondo problemiSiii. We calculate the 
tunneling current through the DQD systems in terms of Green's function techniques for non-interacting systemsiSiiS 
as well as the master equation method. Although electron spin is crucial in the previous theoretical proposals, here we 
disregard it and focus on the quantum interference properties of spinless electrons with/without inter-dot Coulomb 
I 1 interaction. 

This paper is organized as follows. In Sec. ^ a standard tunneling Hamiltonian is employed to describe an AB 
interferometer containing a laterally coupled DQD. We introduce the indirect coupling parameter a. The current 
formula for non- interacting case is provided in Sec. lIIII and the visibility of the AB oscillation is discussed in the large 
bias limit. In Sec. IIVI we provide the current expression in the limit of a strong inter-dot Coulomb interaction. In 
some situations, the current is completely suppressed even when there is a large bias. Our results are summarized in 
Sec. [3 Three sections in the Appendices provides the detailed solutions of the master equation. 
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II. MODEL AND FORMULATION 

We studied laterally coupled double quantum dots (DQD) both of which are tunnel- coupled to left (L) and right 
(R) reservoirs as shown in Fig. (JTJi). The Hamiltonian is TL = Hr + Hdqd + with 

7-Lr = E Y^kctkCuk, (1) 

Hdqd = E £ ( d l d ( ~ t c (d A d B + h.c.) + Ud A d A d B d B , (2) 
<e{A,B} 

H t = E E E [t ( uh<PMkd<+h.c.], (3) 

u£{LR} h (£{A,B} 

where c? vk {c v k) and d^(rf^) represent creation (annihilation) operators of the reservoir v = L/R and the quantum 
dot C = A/B, respectively. We disregarded the spin degree of freedom and we adopt a large limit for the intra-dot 
Coulomb interaction, hence only one level is relevant in each dot. U and t c characterize the inter-dot Coulomb 
interaction and inter-dot tunneling amplitude, respectively. We chose the gauge such that t c is real and positive. 
^l/r X{4'l/r) represents the tunneling amplitude between quantum dot A/B and the mode k in the reservoir L/R. 

The magnetic flux dependence of the tunneling amplitude is t& (<f> v ) /t^f (<f> v ) = exp(=Fi0„), where the upper (lower) 
sign is for v — L(R), and the effective magnetic flux <j) v = 27r$„/ < I>o, which is defined by the flux <f>„ threading 
through the area formed by DQD and the reservoir v and the magnetic flux quantum <i>o = h/e. This Hamiltonian 
also describes the system of a single dot with two relevant energy levels as shown in Fig. QJ>) when t c = and tj> u = 0, 
where U is now interpreted as an intra-dot Coulomb interaction. 

In general, the tunneling current is obtained with the non-equilibrium Green's function (NEC) formalism by 

1 = "kj deTr{(r L -T R )G<(e) 

+ [h(e)T L f R (e)T R }[G r (e) - G a (e)}}, (4) 

where G r (e) and G a (e) are the retarded and advanced Green's function of the DQD, and G < (e) is the lesser Green's 
functioni&i^. The boldface denotes the 2x2 matrix and f v {t) = 1/[1 + e( e ~ p ")/ fc£iT ] is the Fermi distribution 
function where fi u ,ks and T are the chemical potential of the reservoir i>, the Bolzmann constant, and the absolute 
temperature, respectively. The line-width functions are defined as 

r^(e) = 2 7 r^4 ) *(^)4' ) (^)5(e-e^), (5) 

k 

and the off-diagonal component has the following property T v g A = T AB for v = L or R. The wide-band limit 
approximation disregards the energy dependence of T^,. The Green's functions have the following relations: 

G a (e) = [G r (e)}\ G<(e) = G r (e)V<(e)G a (e), (6) 

where S < (e) = i[/i(e)r L + fii(e)T R ] is the self-energy. We have previously discussed the linear conductance for a 
zero offset A = eb — £a = in^S for non-interacting case (U = 0), where the current formula is simpler as follows: 

/ = I J de [J L {e) - f R (e)]Tr{G r (e)r L G a (e)T R }. (7) 

However, if the interaction U is finite, we have to use Eq. Q as demonstrated in the following section. 

By contrast, we also derived a master equation using the method proposed by Gurvitz et alm^, which is appropriate 
for the large applied bias condition: — (J-r. 3> A, T^,, t c , ksT. This approach can only handle the sequential 
tunneling process. We specify the state of a DQD by the occupation number in these dots (N a ,Nb)- The states 
(0, 0), (1, 0), (0, 1), (1, 1) are abbreviated as I — "0", "A", "_B", "D", respectively. p;,;< is the reduced density matrix of 
a DQD after integrating out the reservoir modes in the total system density matrix. The dynamics of the electrons 
passing through a DQD is characterized by the following set of differential equations with an appropriate initial 
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FIG. 1: (a) Schematic diagram of an AB interferometer containing a laterally coupled DQD. The magnetic fluxes threading 
the left and right sub-circuits are <E>l and &r, respectively, and cause the AB effect, (b) Equivalent model of a single dot 
containing two energy levels, which are in the bias window. 
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The functions T^, (e), which describe the tunneling rate of electrons into (out of) dot(s) C, C'> when an electron already 
occupying the DQD, are obtained by replacing in Eq. (JjjJ with e y k — U . In the large limit for the interaction, 
U S> pl — Pr, the tunneling-in process T^-, is absent and we can set pdd = 0. 

The above two approaches are sufficiently general for us to discuss the effect of interaction and interference in the 
transport through a DQD. However, since the line-width function T^, (£ / £'), which controls the strength of the 
coherence of the transport, is strongly dependent on the microscopic model, and we need further simplification to 
grasp the fundamental physics of this system. Here, we define the indirect coherent coupling parameter cv„, which 
was first introduced ini£. The explicit derivation of a v is described in detail inSi. Using ev„, the off-diagonal part of 
the line-width function becomes 



1 AB 



AA^BB 



,=R<? 



(14) 



where the upper (lower) sign is for v — L{R). All the parameters T^-, a v are independent of energy in the wide-band 
limit. We also disregarded the energy dependence of the effective flux induced by changes in the electron trajectory. 
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The parameter cc^or) characterizes the coherent injection into the DQD from the reservoir L (the coherent emission 
from the DQD to the reservoir R). In general a u is a complex parameter but the magnetic flux dependence is 
factorized in exp(=Fz</>„) as shown in Eq. p4|l. \a v \ = 1 corresponds to full coherence and a v = denotes zero 
coherence, corresponding to a situation where the two quantum dots are independently coupled to the reservoir. For 
simplicity, we assume a v is real and positive in the following argument. There has been a detailed analysis of the 
coherence in a metallic reservoir in^. 

Equation (|8I13|I differs from that obtained with a similar method^ and from that obtained with the gradient 
expansion method 24 . Both differ from ours as regards the sign of T^^ with C 7^ C an d the former is missing the 
first term in Eqs. (|12I13() which represents coherent injection of an electron into DQD from the left reservoir. We 
can check that our formula provides a reasonable result in a limiting situation as follows and in the next section. 
Let us consider a symmetric system, namely, T AA — T BB = 7„, zero flux <f> v = 0, and non-interacting [7 = 0. We 
assume complete coherence a R = oll = 1 and therefore T AB = T BA = r y v . We transform from the dot A/B basis to 
the symmetric/antisymmetric (s/a) state basis under the condition of zero offset A = 0. The density matrix is then 
transformed to a new basis as 



Paa = 


\(Pss 


+ Paa 


+ Psa 


+ Pas 


Pbb = 


\(Pss 


+ Paa 


- Psa 


- Pas 


PAB = 


\{Pss 


Paa 


- Psa 


+ Pas 
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\{pss 


Paa 


+ Psa 


- Pas 



with invariant poo & n d Pdd- Then we define state dependent line-width functions: j^/a = 0- a v)lv with +(— ) for 
a symmetric (antisymmetric) state. The master equation for the new basis is 



dp 



(is + la) POO + Is'Pss + la Paa , (15) 

ils + la )Pss + Is POO + la PDD i (16) 

(is + la)Paa + 7a POO + if PDD , (17) 

(if + la)pDD + ifpss + if Paa, (18) 

(if + if + if + if) + 2it c ] Psa , (19) 



dt 
dpss 
dt 

dpaa 

dt 
dp dp 
dt 

dpsa rl 

dt ~ [ 2 

and there is a similar equation for p as . Because of the relaxation term (7^ + "ff + jf + lf~ )/2, the quantum 
coherence term p sa simply disappears from any initial condition for the steady state limit. Therefore, Eqs. (|15I19|) 
correctly describe the independent dynamics of symmetric and antisymmetric channels with state dependent line- 
width functions. It should be noted that when a v = 1, the antisymmetric state is decoupled from the reservoir u, 
j% = 0. This is because of the perfect destructive interference. 

From Eqs. I| 81130 . we obtain the steady state density matrix at t — > 00 by employing the auxiliary relation poo + 
PAA + Pbb + Pdd = 1- Using the result, the steady current is obtained as follows: 

- = ^aaPaa + ^bbPbb + (Taa + ^bb)pdd + Tf B p AB + T§ A p B A- (20) 

e 

III. NONINTERACTING SYSTEM 

First we discuss the system without interaction (U = 0). For simplicity, we restrict to the symmetric coupling 
situation, = 7 and symmetric fluxes 4>L — 4>R = 0/2. The retarded Green's function is 

G r (e) = ( . e ~ £A + 11 _. ± t c + H^e-^+a R e^)\ 
\t c + ^(a L e l 2+a R e e-e B +il J 



5 



For noninteracting conductor, the transmission probability of the electron with energy e is defined as 

T(e) = Tr{G r (e)r i G a (e)r it }, (22) 
which appears in Eq. J7J. The linear conductance G at zero temperature is obtained in the Landauer formula 

G=^T(0), (23) 

where the energy is measured from the (average) chemical potential of the reservoirs. The explicit formula for zero- 
offset, A = 0, is shown in Eq.(18) of Ref« (The definition of the sign of t c is reversed.) The function T(0) for a = 
and a = 1 at zero flux = has following simple physical meaning: 

T(0) = { (e -lp^ + (e 0+ S+^ fOT <* = °- (24) 

1 J2j) for a = 1, V ^ 



where a = corresponds to two independent Breit-Wigner resonances through the symmetric and antisymmetric 
states with line-width 7, while a = 1 represents Breit-Wigner resonance through only the symmetric state with 
doubled line- width 27. It had been shown that the period of AB oscillation of the linear conductance is 2ir when 
t c = and 47r when t c 7^ 0. For non-integer flux, the eo dependence of the conductance shows Fano line shape when 
a = I t. However, this Fano effect is quickly suppressed if a becomes less than 1. 
The current for a finite bias eV = fiL — at T = is 

e f" L 

I = - deT(e). (25) 

'* J PR 



In the limit of a large bias, this can be evaluated by the contour integral and the result for aL — a R = a and A = 
is 

i2+7 2 (l-a 2 sin 2 !) , , 

tt + T 

Now the period of the current oscillation with the flux is 27r independent of t c . At zero flux, <f> — 0, the current is 
e7 independent of a, which is explicitly checked from Eq. H24JI by replacing eo with eo — e and by integrating with e 
for (—00, 00). The current is the sum of e7/2 from symmetric and antisymmetric states for a — 0, and the current is 
e(2j)/2 from symmetric state for a = 1. The energy offset A dependence of the current is shown for various values 
of a in Fig. It should be noted that for sufficiently large offset A ^> 7, the current is independent of a. 

The current for the large bias limit is also derived by the master equation with T^, = T^, . The result is shown in 
Eq. I |A1| I in 1X1 For ait = a L = ot and (p L = (f) R — (f>/2, we have 

A 2 + 4(1 - c? cos 2 IM + 7 2 (1 - a 2 sin 2 f )) 
A 2 + 4(l-a 2 cos 2 |)(i 2 + 7 2 ) 

which provides the same result for A = as that obtained by the NEG method, Eq. I|2()|) . The current is the maximum, 
/ = e7, at (j) — 2nir and the minimum at <f> = (2n + l)ir with an integer n. The visibility of the AB oscillation is 

= 7(0 = 0) - !(<(> = tQ 4 7 2 c 2 

~ I(tf> = 0) A 2 + 4(£ 2 + 7 2 ) ' 1 j 

As shown in Sec. |H] the Hamiltonian also describes the system of a single dot with two levels. Usually, coherent 
injection process to the multiple levels in a quantum dot does not considered. Here we clarify the condition when this 
is justified. By putting t c = and 4> v = in Eq. IjAlfl . we have the current formula: 

'A 2 +7 2 {4- (a L + a R ) 2 } 

The total current deviates from cy, just the sum of the current via each level, ^ej, by the effect of quantum interference 
when 7^ ckr. This effect is maximum if A = and one of the a's is one and the other is zero, where the current 
becomes \erj- This effect of interference vanishes for large offset |A| ^> 7 and the current is ej. This behavior is 
shown in Fig. ??(left). 
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FIG. 2: A dependence of the current with various coherent coupling parameter a = an for left: infinite bias and for right: 
t^L PR ~ 27. t c = eo = 0, oll = and <j> = 0. 



IV. STRONG INTERACTION LIMIT 



Here we consider the case of U — > 00. The general form of the steady current is obtained in [51 For simplicity, we 
restricted ourselves to the symmetrical coupling T V AA = Tg B = j u . In the special case where zero flux <f>„ = 0, we 
obtain from Eq. (|B1|) 

T _ „ 2 7jm ^TI7J + 1 - cr R (30) 



' 27L + A 2 1 27 I ,+(l-a| )7i ? .-2a I ,aR7j: 



When the electron tunneling-out process to the reservoir i? is incoherent, a_R = 0, the current value becomes the 
classical limit, 



' 27^ + 7fl ' 



(31) 



Interestingly, if oil = o;_r, the current has the same value I incoherent as if the coherent transport is absent. In 
both cases, the current value is independent of A and t c . Under general a v conditions, the current value approaches 
^incoherent if the condition A 2 3> 7^+4t 2 . is satisfied. Similar formula is derived for more general situation (asymmetric 
couplings) in Eq. I|B3I) . 

When an — > 1 and A — > 0, the current is completely suppressed even if we supply the system with a large bias. We 
need to keep «l < 1, since an = oil = 1 provides finite current Incoherent- This is evident from the plot of overall 
dependence on an and oll in Fig.|3Ja). This can be understood as the system being trapped in the 'dark-state', which 
in this context means the antisymmetric state that cannot couple to the R reservoir as discussed in the previous 
section. The steady state density matrix in this limit is 

1 1 

Paa = Pbb = ^) Pab = Pba = — ~> Poo = 0, 

which is the density matrix of the pure state: p = \^ a )(^f a \ with the antisymmetric state \^ a ) = ^=(|j4) — \B)). This 

mechanism has been discussed in a triple dot system as a coherent population trapping (CPT) mechanism 2 ^. The 
current in such a system is estimated in [U] and we found the A dependence is similar to Eq. (|3UI) . 

We demonstrate the collapse of current suppression by applying an oscillating electric fiel d 24 ! 26 ' 27 . We evaluated the 
effect of a weak oscillating field A(t) = 5 cos cot for an = 1, and found the leakage current in the lowest perturbation 
in S, 

n = e lR 5H 1 l + Atl+^) 

2(1 - aL ){ 7 ^ + (4t 2 - c 2 ) 2 + 2 7 I(4t 2 +^)Y [6 > 

which is plotted in Fig. [Jb). The current peaks with a value 2 -y R (i-a L ) -y'-^+wt 2 a * a frequency u ~ ±2t c , which 
corresponds to the emission of one photon and the system transits from the antisymmetric state to the symmetric 
state, that allows the electron to leak in the right reservoir. 
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FIG. 3: (a) Left, the a v dependence of the current in the limit of large U and A = 0. The current is strongly suppressed near 
or — 1. (b) Right, photon induced leakage current caused by applying weak oscillating field of frequency ui. The DQD is in 
the current suppressed condition: an = 1, A = 0. 



The flux dependence of the current in the large bias limit is following. We only consider the symmetric configuration: 
1r = 1l = 7 and <\>l = 4>r = <fi/2 and an = 1 and use Eq. (|B1|) . When the offset A is zero, 



e7- 



2^(1 - cos. 



7 2 + (5 - 2a L )t 2 c - ( aL7 2 + (1 + 2a L )t 2 c ) COS0' 
therefore the visibility of AB oscillation of the current is 1. When the inter-dot tunneling t c is zero, 

2A 2 



3A 2 + 2aiA7sin</> + 27 2 (1 — aj, cos< 



where its visibility is 



4a L 7\/7 2 + A 2 



3A 2 + 2 7 2 + 2a L7v /7 2 + A 2 ' 



(33) 



(34) 



(35) 



which is a monotonic decreasing function of | A|/7 and a^. The linear (two-terminal) conductance should be symmetric 
with respect to the reversal of the flux in accordance with Onsager's relation, hence is an even function of <$>. However, 
the current shown Eq. (|34() is obviously not an even function. This is not the problem since we are discussing the 
current in a non-linear regime, that is out of the boundary of Onsager's argument about a linear-response regime. 

We also studied the current expression using the NEG method for the special condition A = 0. For simplicity, 
all tunneling amplitudes are the same and we use the single parameter 7. We transform the basis to the symmet- 
ric/antisymmetric basis with e s = eo—t c and e a = £o + t c and eo = {&A + £_b)/2. The inter-dot interaction Hamiltonian 
becomes Un s n a where n n = (v n d n with n = s/a. On this basis, the line- width function matrix is 



= 7 









1 - OV 



(36) 



We define the total line-width T = T L + T R . We calculated the retarded Green's function Gl^ m (t — tf) = —i0(t — 
t')({d n (t), d' m (£')}), where n, m — s/a using the equation of motion (EOM) approach combined with the lowest order 
decoupling approximation, which is after Fourier transform 



^*nm (^) 



^nmX(e) 



+ 4r„„x(c)' 



(37) 



where x(e) = 1 
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Then we have in the limit of U — > oo, 

GUe) = ^ v {n f, ( ry (38) 

and the lesser Green's function is derived from Eq. |JB}. Finally, the populations of the states are given by 

(%) = / ^G< n (e), (39) 

which should be evaluated self-consistently. The linear conductance at zero temperature for ul — an — a is obtained 
by 

-2(asin(J)) 2 (l - <„.))(! - (n a ))_if^±^_ }; (40 ) 



where 7 s/a = 7(1 ± acos(|))(l - (n„/ s )). 

In the limit of a large bias, Eq. I|39ll reduces to the following self-consistent equation 



and the solutions are 



(l±a L )(l-(n a/g )) 

K/a = 7TT T > ( 41 ) 

' 2 ± a l ± a r 



( n «/o) = o o 1 Y2- ( 42 ) 

3 - 2ctLa R - (a R y 



Putting these in the formula of the current, 

I = E 7 »> ( 43 ) 

n£{s,a} 

7 s/o = e 7 (l±a £ )(lTafi) 1 ~ (n ° /s) , (44) 
we finally obtain the same result as that obtained by the master equation, Eq. (|30|l with A = 0: 

I = 2e 7 2 . (45) 

o — 2aL«i? — cti? 



V. CONCLUSIONS 



We discussed the effect of quantum interference on the transport through a quantum dot system. We stressed 
the role of the indirect coherent coupling parameter a, which provides constructive/destructive interference in the 
transport current depending on its phase. We derived the current using the non-equilibrium Green's function method 
as well as the master equation method in the sequential tunneling regime. For a large inter-dot Coulomb interaction, 
the current is strongly suppressed by the quantum interference effect, where the current is restored by applying 
oscillating resonant field. 
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APPENDIX A: DERIVATION OF CURRENT FOR NON-INTERACTING SYSTEM 

Steady state density matrix elements are derived from the master equation by setting = in Eqs. 1(8113(1 . When 
the interaction is absent, T^, = r^-,. Since the six algebraic equations are not independent, we need the equation 
for conservation of probability, poo + Paa + Pbb + Pdd = 1- Evaluating the current with Eq. 1(21) (1 . we obtain 

Nj 

I = e 7 ^, (Al) 

Nj = A 2 + 7 2 [4-2(a 2 L +a 2 B ) + a 2 L a 2 R sm 2 (0 L +(f> R )} 

+ 2i 2 [2 — a 2 L cos 2 4>l — a R cos 2 cj) R ], 

Di = A 2 + 7 2 [4 - a 2 L - a\ - 2a L a R cos{<t> L + 4>r)] 

+ t 2 [4 - (ctL cos c/)l + or cosier) 2 ]. 



APPENDIX B: DERIVATION OF CURRENT IN THE STRONG INTER-DOT INTERACTION LIMIT 

We solve the master equation Eqs. I|81 13(1 for U — > oo, namely, neglecting the term pdd and T^, in the steady state 
condition. First, we restrict ourselves to the symmetrical coupling T AA — T BB = 7 „ and the result is 

I = e^g^{A 2 ~ 2a 2 R t 2 c cos2<^ + (1 - a 2 fl ) 7 2 + 2(2 - a 2 R )t 2 c }, (Bl) 

D = A 2 {2 7L + 7R ) + 2a L a R Aj L j R sin((t) L + (f) R ) + (l~a 2 R )j R 

-2a L a R j L {j R + 2t 2 c ) cos(</> L + (j) R ) - Aa L a R ^ L t 2 c cos(</> fl - cf) L ) 

-2a 2 RlR t 2 c cos2^ + 2(2 - a 2 ) lR t 2 c + 2 7l ( 7 2 + \t 2 L ). (B2) 

In the model of a single quantum dot with two levels (t c = 0,0,, = 0) as shown in Fig. db), current is evaluated 
for the most general choices of T^, , 

I -^incoherent 

K(A), (B3) 

(^AA +^Bb)^AA^BB j-g^-j 



^incoherent 



pfl p-R 1 p-R pL 1 pfl pL : 
1 AA BB ' AA BB ' BB L AA 



K(A) = 



4A 2 + (i- a 2 f ,)(r^ + rg s ) 2 



1A ; 1 {±^ME2A£BB±ElA£BBj£BBElAz3s^^ (rR +r« V 

r*f TrWZTrx .rk^4-rk^r L . . \ l AA' l bbj 



where n = y r^T^r^r^. In this model, the incoherent current is suppressed if one of the coupling to the right 

reservoir, T AA or T BB , is very smallS^. The corresponding local state in the quantum dot is now the 'dark state' to 
suppress the current. The function K(A) is 1 when the coherent coupling in the right reservoir is absent, a R = 
irrespective of the values of A, ol, and so forth. When a R — » 1, K{A) is suppressed for small A, while for large offset, 
A, K{A) 1. 



APPENDIX C: DERIVATION OF CURRENT THROUGH THE TRIPLE DOT SYSTEM 

We consider the triple-dot model used in Refi^ with two lateral dots, A and B, with energies ea and e R coupled to 
left reservoirs independently and one dot, C, with energy e c coupled to right reservoir. The dot A and B are tunnel 
coupled to C with amplitude t. Because of the charging effect, the total number of electrons is zero or one and the 
applied bias is very large. Setting up the master equation as done in the main text, we obtain the formula of steady 
current: 

2 7L7fl A 2 



{2 7L + 7R + 7>(2(e - £c) 2 + ^±%A2 + 8 7L t 2 ' 
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where eo = £A + £g and A = £g — ea- This formula resembles the result of Eq. (|30l) with — 1, where the current is 
strongly suppressed near A = 0, while the behavior for large |A| is different. 
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